Entanglement Entropy of the Early Universe in Generalized Chaplygin Gas Model 
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We provide an explicit calculation of the evolution of the cosmic entanglement entropy in the 
early universe before the matter dominant era. This is made possible by invoking the generalized 
Chaplygin gas (GCG) model, which has the advantage of preserving unitarity and providing a 
smooth transition between the inflation epoch and the radiation dominant era. The dynamics of 
the universe is described by the quantization in the minisuperspace of the GCG model, following 
the prescription proposed by Wheeler and DeWitt. Two sources of contribution to the cosmic 
entanglement entropy are considered: one from the homogeneous background where the observable 
and the unobservable regions of the universe are entangled and the other from the inhomogeneous 
cosmological perturbations where different modes are entangled. We find that the homogeneous 
contribution grows exponentially at the very beginning of the inflation, but decreases during the 
radiation dominant era. Conversely, that from the cosmological perturbation is found to decrease at 
first and then increase after reaching a minimum value. The net result is that the total entanglement 
entropy reaches a minimum at an early stage of the inflation and then increases throughout most 
of the inflation and the entire radiation dominant era. 

PACS numbers: 98.80.Bp, 98.80.Qc, 89.70.Cf, ll.15.Bt 
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I. INTRODUCTION 



An increasing amount of observational evidence sup- 
ports the notion that the early universe has undergone 
an epoch of inflation. However we are still far from under- 
standing the underlying assumptions and resolving some 
of its most crucial issues. One of the open questions 
about inflation is the entropy problcm[l]^]^3j], that is, 
why did the universe start from an extremely low en- 
tropy initial state, which in turn gave rise to the arrow of 
time? Since the initial entropy evaluates the probability 
of spontaneous formation of a homogeneous domain of 
an inflationary universe, the very low initial entropy im- 
plies the unnaturally small probability for our universe 
to evolve to the current state through inflation. Tanta- 
mount work has shed lights on possible solutions to this 
problem. For example one approach is the spontaneous 
eternal inflation where inflation can occur both forward 
and backward in timeQ. This has also been treated in 
string theory in the form of bubble cosmology ^4] where 
universes that realize the complete set of string vacuum 
solutions were constructed ^SjQ. More recently the tri- 
partite partition of a physical systemQ was introduced 
where the entropy was found to decrease dramatically 
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after inflation as a means to solve the problem. 

We note, however, that the definitions of the entropy 
budget are hugely different in these various approaches. 
For example at late times when local thermal equilibrium 
can be reached in the universe, the notion of thermal en- 
tropy is applicable in issues related to cosmological struc- 
tures such as galaxies or clusters. For strongly gravitat- 
ing systems such as black holes, holographic entropy 
satisfying the Beckcnstcin bound is used and most works 
approximate the total entropy of the universe as the sum 
of that of all black holes Q IS 11|; in string theory Q 
the Weyl curvature hypothesis [l| is often assumed which 
suggests that the arrow of time is due to the vanishing of 
the Weyl tensor that effectively acts as the entropy; for a 
universe dominated by ultrarelativistic matter, the total 
entropy coincides with the total number of particles in 
the universe _12J. Apart from this lack of conformity in 
the definition of entropy when solving the entropy prob- 
lem, all previous work only evaluated the entropy qual- 
itatively due largely to the lack of a working model to 
track the entropy evolution in the early universe. 

In principle the thermal entropy in a purely gravitat- 
ing system is ill-defined for lack of a universal tempera- 
ture. Also, statistical mechanics is not directly applicable 
due to the long-range nature of gravitational force. Self- 
gravitating systems thus possess unusual features such as 
the absence of the global entropy maxima, different from 
conventional thermodynamical systems[13j. Even so, the 
^f^ss of information in a quantum gravitating system can 
still be evaluated by von Neumann entropy, which mea- 
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sures the degree of bi-partite entanglement. Since the 
early universe can be described as a quantum system 
without strong decoherence thereafter, it is justifiable to 
invoke von Neumann entropy in finding the evolution of 
the total entropy of the universe, which is a dynamical 
quantum system far from thermal equilibrium. We will 
thus use the entanglement entropy instead of entropy in 
the following discussion to help clarify our setup and to 
quantify the amount of information in a quantum system. 

In order to track the entanglement entropy of the uni- 
verse, a specific model is needed. Most existing theories 
about the early universe necessitate the reheating process 
to connect the inflation to the subsequent radiation dom- 
inant erafl^. However typical reheating mechanism such 
as the decay of the scalar field is not adequate to render 
tangible calculations of the entanglement entropy. As a 
remedy, we resort to an alternative approach, i.e., the 
generalized Chaplygin gas model (GCG), which bridges 
the inflation epoch and the radiation dominant era via 
a single equation of state and therefore a smooth tran- 
sition. Though a phenomenological theory, GCG repro- 
duces well the CMB power spectrum[l^. Although we 
invoke GCG so as to accomplish tangible calculations on 
entropy, the method developed in this work is general and 
can be used in arbitrary set up. The underlying effective 
Lagrangian of the theory, which is consistent through- 
out the whole period, yields a solvable Wheeler-DeWitt 
equation and the background evolution, both of which 
are crucial for the evaluation of entanglement entropy. 

In Sectini we briefly introduced the GCG inflation 
model used to approximate the evolution of the early 
universe. This is the best candidate model so far in re- 
gards to providing a manageable framework for a speciflc 
calculation and the ability to reconstruct the existing ob- 
servational data and make important predictions [ISJ. We 
expand the detailed method of finding the entanglement 
entropy in Sect lllll and llVl which make up two parts of the 
whole entanglement entropy: the homogeneous part and 
the inhomogeneous part. Numerical results are presented 
as figures followed with brief analysis in discussion where 
a decrease in homogeneous contribution is seen in radi- 
ation dominated era and a saturation to a lower bound 
appears in perturbation entanglement entropy. 



II. GENERALIZED CHAPLYGIN GAS 

When applied to cosmology, the Chaplygin gasp^ 
models the change of cosmic content by regulating the 
equation of state of the background fluid instead of 
the form of the potential. It was flrst suggested by 
Kamenshchik 16] in an attempt to smoothly interpolate 
the de Sitter phase and the radiation dominant era with- 
out ad hoc assumptions. Generalized Chaplygin gas 
models (GCG) were subsequently suggested[17]U8],19], 
some of which managed to unify dark energy and dark 
matter [20II. GCG was also introduced, within the frame- 



work of FRW cosmology, as an exotic background fluid to 
admit supersymmetric generalization [23|- The equation 
of state of GCG is given bvflU: 



A B 



(1) 



From the conservation of energy the pressure for the 
Chaplygin gas is 



1 + ^-4(1 + a)/ B 
3(1 + a) 



One can therefore flnd the underlying minimal coupling 
scaler field (/> and its potential 1^(0) through the relations 
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We require that the inflation takes place before radiation 
dominant era, which leads to the restriction: 
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4(l + a) l + P' 



(4) 



Two more conditions should be imposed on a and /? to 
have a model that interpolates between an early inflation- 
ary phase of the type of quintessence and a subsequent 
radiation dominant phase: (i) the energy density must in- 
duce a period of inflation and (ii) the inflation should not 
cause a super-inflationary expansion; i.e. < H. That 
is, a super-accelerating phase of the universe, where the 
energy density grows as the universe expands, should not 
occur. By combining these two ansatz with the above in- 
equality (Eqdl])), we can easily deduce that the set of 
allowed values of a and /3 satisfy 
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1 -I- a < 0, 
1 + I3<0, 
a) <1 + p. 



(5) 



The statistical feature of the Chaplygin gas can be de- 
duced from the flrst law of thermodynamic at much later 
time when thermal temperature can be well-defined: 



TdS ^ d{p{T)V) + p{T)dV, 
and its equation of state. 



The entropy is then 



p = wp. 
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Let the holographic entropy be Sa — nr\/G. The evo- 
lution of the total entropy can then be parameterized as 
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Since the early universe can be treated as a quantum sys- 
tem without significant decoherence afterwards, thermal 
entropy should be replaced by the von Neumann entropy, 
which measures the degree of quantum entanglement. 
Without relying on further assumption of the area-law 
behavior of the holographic entanglement entropy, we try 
to find the total entropy budget by resorting purely to 
the von Neumann definition. 

Once a quantum system is entangled, the wave- 
function of the system under consideration cannot be 
factorized. From the observational point of view, this 
means that we fail to access the entire phase space of 
the system since any single measurement will destroy the 
entangled state into a factorized state so as to conceal 
its entangled counterpart, which is exactly the process 
of decoherence. As a result, the degree of entanglement 
measures the loss of information, which is equivalent to 
the von Neumann entropy in a quantum system. 

As we know, entanglement entropy from the von Neu- 
mann measurement counts the degree of bipartite entan- 
glement and the result is dependent on the division of the 
system. Even for a system which is pure, the entangle- 
ment entropy of one subsystem can measure its degree 
of entanglement with the other subsystems. Multipar- 
tite entanglement entropy is thus needed to evaluate the 
overall degree of entanglement between all possible pairs 
of subsystems. There are definitions of multipartite en- 
tanglement entropy such as the global entanglement that 
maximizes the fidelity between a quantum state and a 
separable state[2^. This definition, however, is difficult 
to employ for general systems. 

We suggest that one reasonable way to quantify the 
total entanglement entropy of a quantum system is to 
sum up von Neumann entanglement entropies of all sub- 
systems. A common approach to estimating the total 
entropy of the late time universe [9], for example, is by 
adding the entropy of each black hole, which follows this 
same principle. The only difference is that this approach 
only considers the entanglement entropy between an in- 
dividual black hole and the rest of the universe, whereas 
we here propose that the entanglement entropy between 
every pair of conceivable subsystems should be accounted 
in order to find the total entanglement of the whole quan- 
tum system. The value of multi-partite entanglement en- 
tropy in our definition can be infinitely large. However 
the absolute value of the total entanglement entropy may 
be less essential than its change. By focusing solely on 
the change of the entanglement entropy, our method can 
serve to unveil the dynamical behavior of the entangle- 
ment of the quantum system. One other virtue of our 
generalization lies in that the important features of the 



entanglement entropy found in previous studies such as 
subadditivity are well preserved. 

Let us first define the early universe as a quantum sys- 
tem. The inflation in the early universe can be roughly 
divided into four stages: the onset of the inflation as a 
quantum gravitating system, the scaler fleld perturba- 
tions, the horizon exit and the radiation dominant era. 
Since the horizon re-entry happens after the radiation 
domination, the perturbations evolve during the radia- 
tion dominant era purely due to the expansion of the 
background metric. Since entanglement is not a statisti- 
cal quantity, the entanglement entropy can be evaluated 
in systems far from equilibrium, which is the case for the 
highly dynamic universe during inflation. The evolution 
of the entanglement entropy in these different stages of 
inflation will be calculated explicitly in the following sec- 
tions. 

Suppose we can describe the state of a system com- 
pletely by a set of discrete variables O = Oi,02,---- 
Once we find the variables of the system in certain values 
Oi, the system is in a pure state so that the conform- 
ing entanglement entropy is zero. However, if we know 
only the probability distribution Pq, then following von 
Neumann's definition, the entanglement entropy becomes 
S = —^^PolnPo- If the variable is continuous, one 
takes the continuum limit and the entanglement entropy 
is in the form S = -J P{0) \nP{0)dOQ^. 

In order to evaluate the entanglement entropy based 
on this generalized multi-partite formalism, we need to 
determine the observables and their probability distri- 
butions. During the inflation epoch, it is reasonable to 
use the cosmic scale factor a and the perturbation ampli- 
tude 5k as our observables, since other physical quantities 
can either be determined by these two or are unobserv- 
able under current consideration. And there are two dis- 
tinct contributions to the entanglement entropy of the 
early universe. We find the entanglement entropy con- 
tributed by a in Section lHll bv solving the wave-function 
from Wheeler-DeWitt (WDW) equation and then inte- 
grate it over a. This part of the entanglement entropy is 
due to the evolution of background metric and is thus a 
homogeneous contribution. On the other hand, the inho- 
mogeneous part, accounted in Scction llVl is contributed 
from the entanglement between two opposite modes k 
and —A: of the adiabatic perturbation in a single-field 
Chaplygin gas. The total inhomogeneous entanglement 
entropy is then obtained by integrating over the entire 
mode spectrum. 
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III. ENTANGLEMENT FROM THE 
WAVE-FUNCTION OF THE EARLY UNIVERSE 

Based on the canonical quantization approach to quan- 
tum gravity, WDW equation iJ^* = o[2^[2^ determines 
the geometry of space-time of the minisuperspace of a 
FRW universe fihed with Chaplygin gas. For a closed 
FRW model the minisuperspace Lagrangian is[27j(28j[29j 



L = ^y/^R + LcG- 



(10) 



The early universe is assumed to be homogeneous and 
isotropic that can be well approximated with the metric 



ds2 = -N^dt^ + a{tf 



dr2 



l-r2 

We can then obtain the Lagrangian as 



(11) 



Clearly, the wave-function of the universe depends on 
the canonical variables a and r: \E'(a, r), where r is the 
conformal time. Following the standard procedure, we 
solve the wave- function by replacing Pa with in the 
Hamiltonian, and the WDW equation becomes, 



12 da^ 



-3a' + —{ — 



A 



B 



^{a,T) = 0. 
(17) 



The general solution to the wave-function can be ex- 
pressed as 



(18) 



where 



^CG 



1 



-gR 



A 



N 
B 



Na , (12) 



1/(1+q) 



(13) 



^^ = eM-S\a)]. (19) 
Ci and C2 can be fixed by the boundary condition, and 
S'*(a) satisfies the Hamilton- Jacobi equation: 



The momentum conjugate to the scale factor a, by defi- 
nition, is 



Pa = 



d{LcG 



-9R) 



da 



6aa 



(14) 



The Hamiltonian of the homogeneous and isotropic FRW 
universe with a positive curvature is thus 



(15) 



In our parameterization, the effective potential is 
-4 / A B 



l/(l+«) 

I 



da 



= 121/0 (a) 



(20) 



The exact analytic solution for the effective potential in 
the form of Eq ([T6|) is difficult to find. However, at very 
early times where the scale factor is very small, the po- 
tential can be approximated by 



Via) = 1 



(3 + 1 \ _M± 



6(a + l) 



a °+i . 



(21) 



Then the wave-function has a relatively simple solution: 



^ = exp 



-2- 



^/2a(^)(a-H)V ^^^^"'fff°^*^- 



a + 1 



V2ao(c 



1 



r(6a-;3+5)a( 



3+1 

■ Q + l 



a+1 



2a- (3 + 1 



2a -(3 + 1 



(22) 



We take uq — 10 as our initial condition. The entropy 
can then be obtained by tracing over the even horizon: 

S = Tr(plnp) (23) 



a(r) 



1'**ln(**«')da + S'o, 



(24) 



where the wave-function is already normalized and Sq 
refers to the initial entropy. In this case, time variable is 



replaced by the scale factor and the background metric is 
determined by the equation of motion of the Chaplygin 
gas. 



Carrying out the integration, the entropy is found to 



be 
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S = fia,/3,B,A)T 



4a- (3 + 3 
2a- P + r 



a{ty 



1) 



2a- (3 + 1 



- (6a-;3+5) 



Q+1 



Sid, 



(25) 



v 



where the coefficient f{a,(i,B,A) is determined by 
a, (3, B and A exphcitly. The numerical result is shown in 
Fig. [T] We see that the entanglement entropy increases 
exponentially at the first stage of the cosmic inflation 
when the scale factor is also under exponential growth. 
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FIG. 1: The evolution of the entanglement entropy con- 
tributed from the entanglement between the observable uni- 
verse and that outside the horizon froma = 10~® to a = 10^, 
obtained by solving the WDW wave-function. 



At a later time during the radiation dominant era, the 
effective potential can be approximated as 



V{a) - -B'/^^+-'> 
^ ' 3 



(26) 



and the wave-function is again solvable, and we find 



exp 











ait)) 









(27) 



The corresponding entropy is found to be 



S 



3ao 



exp 



2B2(<>+i) 



V3i 



{2B^^ - \/3aoSWTT^ Ei 



2B2(°+i) 



\/ia{t) 



(28) 



-\/3B2(-+i) floexp 



2i32(=+i) 
\/3ao 



\/3aoS5(^ - 2B^\ Ei 



a{t) exp 



2B2(c+i) 



2B2(c+i) 
^/?,a(t) 



\/3ao 



This solution can be expressed in a simpler form: 



S = -ifiEi 



2B2(<>+i) 
y/ia{t) 



K2a(t) exp 



2B2(<>+i) 



(29) 



where Ei(x) stands for the exponential integral function 
which increases monotonically for a; > 0, and ifi, K2 and 
are positive constants. Interestingly, we find that the 
first derivative of the entanglement entropy is negative 



and thus it decreases monotonically. This conclusion is 
verified by a numerical computation (see Fig. [2]) 



However, since the total entanglement entropy con- 
tains other components apart from this homogeneous 
contribution from the wave-function of the background 
evolution, such as that from the cosmological perturba- 
tions, this decrease in entropy does not mean that the 
total entropy of the universe necessarily decreases. 
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The field iy{k) satisfies in the Fourier space the equation 
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FIG. 2: The evolution of the entanglement entropy con- 
tributed from the entanglement between the observable uni- 
verse and that outside the horizon in the radiation dominant 
era from a = 10^^ to a = 10^^. 



IV. ENTANGLEMENT ENTROPY OF 
COSMOLOGICAL PERTURBATION 

The Mukhanov-Sasaki variable v represents the linear 
combination of the inflaton perturbations of the gravita- 
tional potential. The comoving curvature perturbation 
is then Ci^jt) — v{x,t) / z{x,t), with z{x^t) — a^/AG 
and e = —dlnH/dlna. The conjugate momentum of C 
is denoted as tt = dp,C- The Gaussian random state in 
a single scalar field inflation is then characterized by the 
covariant matrix between two modes k and -k, which is 
related to the density matrix asfsoj 



C = -TripV,V^) = { 



V ■ 



Each component can be identified as 

3 

The corresponding entropy can then be found: 
S = 2[{n + l) ln(n + 1) - n In n)], 



(30) 

(31) 
(32) 

(33) 

(34) 
(35) 



V —V ^ 



-V = 0, 



(36) 



which governs the evolution of C,{h). 

So the basic steps to calculate the entropy are: i) evolve 
the background metric dictated by Eq. ((38|) : ii) solve the 
comoving curvature perturbation at time t for different k 
according to Eq. p6p ; iii) integrate the total entropy over 
the momentum space following Eq. p9p . The numerical 
results is shown in Fig. [3] 
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FIG. 3: Numerical result for the entanglement entropy con- 
tributed from the cosmological perturbation from the initial 
state of asymptotic Minkowski perturbation at a = 10"^" to 
a time well after the horizon exit for most modes at a = 10"^*. 



Mind that though from the current observation, we 
cannot determine the exact starting time for the expected 
cosmological perturbation, a rough upper limit can be 
given by requiring the energy scale not to exceed the 
planck mass. The corresponding time for this limit is 
at approximately a = 10""^^ much below our concerned 
value in calculation. 

One way to look into the entanglement behavior at 
much earlier time is by assuming the vacuum state as the 
asymptotic initial condition. Here, for example, the free 
vacuum is heuristically taken to be the Bunch-Davies vac- 
uum. Inserting the background evolution into our GCG 
mo del [ill, we find 



2(1- 



4 2(l + „) . / — t 

a) V 3 



1+^ 



(37) 



where k} = SttG, so that the entropy of the whole space 
under the assumption of Bounch-Davies vacuum is 
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dq ■ 2 



ic; 



ic; 



In 



IC 



t^-'q^ 



(38) 



t^-'^q^ 



where ^max = ^ir/lp and the particle horizon of the sys- 
tem, L = j^dt/a{t), determines the ultraviolet cutoff 
gmin = 27r/L, which can also be expressed in terms of t 
in our model as Qmin — [27r7''/(l — v)\t. 

It is easy to prove that for 
\\5l\H'^-''ql^J{'i'KGv^-1^'')\ < 1/2 + 1/e, the entropy 
increases monotonically with time. Since the scalar 
perturbation reaches its asymptotic maximum at the end 
of the free scale perturbation, i.e., |(5inax| = | <^horizonoxit | , 
we can give a rough estimate of the time when entropy 
henceforth increases incessantly: 



t > 



(l + 2e-i)Gz/3(l-jy) 



l/3(i^+l) 



(39) 
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approach, we have evaluated the entropy evolution con- 
tributed from the homogeneous background and the inho- 
mogeneous perturbation dynamics, from the very begin- 
ning of the inflation to the radiation dominant era. We 
found that the entanglement entropy contributed from 
the wave-function in the radiation dominant era under- 
goes a monotonic decrease. But the total entanglement 
entropy still increases during the same period due to the 
larger counteracting contribution from the curvature per- 
turbation. However we found, as shown in FiglU that at 
the much earlier time when the inflation has just started, 
the total entanglement entropy has indeed reached a min- 
imum at one time. This interesting turnover happens at 



1 in our unit, which corresponds to a 
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DISCUSSION 



the usual notation. This unexpected result underlies the 
possible existence of a lower bound in the entanglement 
entropy that conforms with the proven bound in a much 
simpler quantum system[3^. Whether the existence of 
such a lower bound is universal remains to be tested. 



We have seen that the generalized Chaplygin gas 
model, though a phenomenological approach, may well 
provide a smooth transition between the inflation and 
the radiation dominant phases in the early universe. On 
the other hand, it fails to capture one essential feature 
common to most other inflation models, namely the post- 
inflation reheating process where copious particles are 
created. However, if as suggested ^3] [32] that the evolu- 
tion of the universe is unitary, GCG model which pre- 
serves unitarity may have the advantage over other infla- 
tion models in solving the entropy problem. 

One should keep in mind that the results obtained in 
the above sections, in particular that shown in Figll] 
m and m only concern the variation instead of the ab- 
solute value of the total entanglement entropy. In our 
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